Generic features of modulational instability in nonlocal Kerr media. 
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The modulational instability (MI) of plane waves in nonlocal Kerr media is studied for a general, 
localized, response function. It is shown that there always exists a finite number of well-separated 
MI gain bands, with each of them characterised by a unique maximal growth rate. This is a general 
property and is demonstrated here for the Gaussian, exponential, and rectangular response func- 
tions. In case of a focusing nonlinearity it is shown that although the nonlocality tends to suppress 
MI, it can never remove it completely, irrespectively of the particular shape of the response func- 
tion. For a defocusing nonlinearity the stability properties depend sensitively on the profile of the 
response function. It is shown that plane waves are always stable for response functions with a 
positive-definite spectrum, such as Gaussians and exponentials. On the other hand, response func- 
tions whose spectra change sign (e.g., rectangular) will lead to MI in the high wavenumber regime, 
provided the typical length scale of the response function exceeds a certain threshold. Finally, we 
address the case of generalized multi-component response functions consisting of a weighted sum of 
N response functions with known properties . 
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1. Introduction 



The phenomena of modulational instability (MI) of plane 
waves has been identifiecLand studied in various phjssi 
cal systems, such as fluidslll, plasmao, nonlinear opticaJ'l 
discrete nonlinear systems, such as molecular chain; 
and Fermi-resonant interfaces and waveguide array; 
dispersive nonlinear directional couplers with the change 
of cefractive index following a exponential relaxation 
lawLj etc. It has been shown that MI is strongly af- 
fected by various mechanisms present in nonlinear sys- 
tems, such as higher order dispersive terms in the case 
of optical pulsesa, saturation of the p^nlinearityO, and 
coherence properties of optical beamsEEl. 

In this work we study the MI of plane waves propagat- 
ing in a nonlinear Kerr type medium with a nonlinearity 
(the refractive index change, in nonlinear optics) that is 
a nonlocal function of the incident field. We consider a 
phcnomcnological model 



id z tp + -ypti> + sip 



R(x' - x)\ip(x')\ 2 dx' = 0, 



(1) 



for the wave modulation ip where x is the transverse 
spatial coordinate and s = 1 (s = — 1) corresponds to a 
focusing (defocusing) nonlinearity. The evolution coor- 
dinate z can be time, as for Bose-Einstein Condensates 
(BEC's), or the propagation coordinate, as for optical 
beams. We consider only symmetric spatial response 
functions that are positive definite and (without loss of 
generality) obey the normalization condition 



R(x)dx = 1. 



(2) 



Thus we exclude asymmetric effects, such as those gener- 
ated by asymmetric temporal response functions (with x 
being .time), as in the case of the Raman effect on optical 
pulsest! 

In nonlinear optics (|lj) represents a general phe- 
nomenological model for media in which the nonlinear 
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refractive index change (or polarization) induced by an 
optical beam is determined by some kind of a transport 
process. It may include, e.g., heat—conduction in ma- 
terials with a thermal nonlinearityEjjLj or diffusion of 
molecules or atoms accompanying nonlinear light prop- 
agation in atomic vapourstLa. Nonlocality ralsct, accom- 
panies the propagation of waves in plasmat£rE3, and a 
nonlocal response in the form ([!]) appears naturally as a 
result of many body interaction processes in the descrip- 
tion of Bose-Einstein condensatesEil. 

The width of the response function R(x) relative to 
the width of the intensity profile \ip\ 2 (x,z) determines 
the degree of nonlocality. In the limit of a singular 
response we get the well-known nonlinear Schrodinger 
(NLS) equation, which appears in all areas of physics. 
Here the focusing case (s=l) produces MI of the finite 
bandwidth type, while the defocusing case (s=— 1) pre- 
dicts modulational stability!! When the width of the 
response function is finite but small compared to the in- 
tensity, the modei, ([j]) is approximated by the modified 
NLS equationE3-Ea 

id„il> + \d 2 x ^ + sU\ 2 + 7 ^|vf]V> - 0. (3) 
Here 7 <C 1 is defined as the second virial of R(x) as 

^ poo 



1 ; 2 



(6) 



7 



x R (x) dx 



(4) 



and it scales as 7 ~ a 2 , where a is the characteristic 
length of the response function. In contrast to the local 
NLS limit (7=0), the MI now depends not only on the 
sign of s, but also on the intensity of the plane wavesEH. 
Finally, in the case of strong nonlocality it has been 
shown that (^) simplifies to, a linear model, and hence 
there is no MI in this limitEZl. 

MI has thus been studied in different limits. The gen- 
eral case (|l|) has recently been investigated with respect 
to MI and compared with the weakly nonlocal limited. 
Here we present an analytical study of the full nonlocal 
case with arbitrary profile R(x) whose spectrum obeys 
a sufficient degree of smoothness, with particular em- 
phasis on generic features of the MI. The present-paper 
complements and extends the results obtained ii£3. 



2. MI in the nonlocal NLS equation 

The model (Q) permits plane wave solutions of the form 

ip(x, z) = yfp^ exp(ik x - iuj z), p > 0, (5) 

where po, fco, and luq are linked through the nonlinear 
dispersion relation 



Following!!, we perturb the plane wave solutions (||) - 
(^|) as follows: Assume that 

ip (x, z) = + u (£, r) + iv (£, t)] exp(i6> ) (7) 
£ = x — koz, t = z, (8) 
6»o = k x - uj z, (9) 

where u and v denote the real- and imaginary part of the 
perturbation. Inserting this expression into the nonlocal 
NLS equation (Q) and linearizing around the solution (||) 
- (||) yields the equations 

d T u + ^d' 2 v = 0, 

d T v - \d 2 u - 2 SP0 Ho R(C - 0«(£'> r)< = 0, 

(10) 

for the perturbations u and v. By introducing the 
Fourier transforms 



u(k, r) 
v(k, t) 
R(k) 



u(£,r)exp[jfc£]ci£, 
v(£,t) exp d£, 
i?(£)exp[z^]d£, 



(11) 
(12) 
(13) 



and exploiting the convolution theorem for Fourier trans- 
forms, the linearized system is converted to a set of or- 
dinary differential equations in /c-space 



d r u — ^k 2 v = 0, 

d T v + \k 2 u — 2spoRu = 0, 
which can be written in the compact matrix form 

d T X = A2L, 
where the vector X and matrix A arc defined as 



x = 


u 


> 4 = 




V 









\k 2 



■\k 2 



2sp R{k) 
The eigenvalues A of the matrix A are given by 



A 2 = k 2 Po 



sR(k) k 2 

4po . 



(14) 



(15) 



(16) 



(17) 



The general dispersion relation (|T^) constitutes the basis 
of our study of MI. 

Let us summarize the properties of the spectrum R(k): 
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1. Since R(x) is realvalued and symmetric, then so is 
R(k), i.e. R(k) = R(-k) = R*(k). 

2. The normalization (|^) implies i?(0)=l. 

3. The symmetry condition for the Fourier transforms 
imposes R'(Q)=0, i.e. the spectrum has a critical 
point at k=0. Here and in the following prime de- 
notes differentiation with respect to the argument. 

4. The normalization (||) for R(x) means by as- 
sumption that R(x) is absolute integrable, and 
hence by Riemann - Lebesque lemma we have 
lim,.^ R(k) = OB. 

5. The functions R(x), xR(x) and x 2 R(x) are as=. 
sumed to be absolute integrable, and thus byEj 
R(k), R'(k), and R"(k) are continuous for all k. 

6. The response functions are characterized by typi- 
cal widths and scaling lengths a and they assume 
the generic form R(x) = a~ 1 ^(x/a), where $ is a 
non-dimensional scaling function. 

The spectrum R (k) can be expressed in terms of the 
Fourier - transform $ of the scaling function <3? as 



R{k) = $(crfc) 



$ (C) exp [iakQ d(. (11 



Notice that the list of general properties 1-6 of the spec- 
trum R carries over to $. 

The dispersion relation (|l7]) is now conveniently 
rewritten as 



A 2 /Aq = k 2 (j} (k; s, a) , 
by means of (Eq), the definition 



A = 2p 
the redefinitions 

2\/A)cr -> (j, 



and 



4>(k;s,a) = s$(ak) 



(19) 
(20) 

(21) 
(22) 



The parameter a which measures the width of the re- 
sponse function, plays the role as a control parameter. 
Observe that $ is a continuous differentiable function of 
k. The crucial point in the stability analysis is the prop- 
erties of the function <p. More precisely, by appealing to 
the list of properties 1-6 we can characterize the set Vl of 
k fulfilling the inequality (f>(k;s,a) > for a given value 
of a as follows: 



• If <fi (k) < for all k > 0, then Ct is empty. In this 
case we have modulational stability. 

• If the transversality condition (j)' ^ Ois satisfied at 
the zeros ofcj), the number of zeros is finite. In addi- 
tion, these points are distinct and isolated. In this 
case n is given as a finite union of well-separated, 
closed, bounded subintervals of the positive k - axis. 
In this case we have MI of the finite bandwidth 
type. 

• The breakdown of the transversality condition for 
certain values of a, i.e., 



4> (k) = (j)' (k) = for k — k c , a = a c 



(23) 



describes bifurcation phenomena like excitation, 
vanishing, coalescence, and separation of MI 
bands. A careful analysis of the behavior of <fi in 
the neighborhood of the bifurcation point will re- 
veal the type of phenomena which takes place. The 
number of zeros of <j) will change as a passes the 
critical value a c . The case </>" (k c ) > for a = a ci 
can be interpreted as a condition for two MI bands 
to merge together or that one single MI band sep- 
arates into two gain bands, while 4>" (k c ) < for 
a = a c represents excitation or vanishing of an MI 
band. The critical value a c and the type of bi- 
furcation can be determined as follows: Solve the 
equation 



z$' (z) =2$(z), z>0 



Then the widths a c can be expressed as 



(24) 



s$ (z) 



(25) 



provided s$ (z) > 0. The second derivative of cj) 
evaluated at the bifurcation point can also be ex- 
pressed in terms of z : 



z 2 $" (z) 



(26) 
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From now on we will assume that the transversality 
condition is satisfied at the zeros of <j), except for a set 
of isolated bifurcation points. This means that the be- 
havior of the graph of spectrum $ (ak) relative to the 
parabola fc 2 determines the stability properties. We will 
discuss this aspect in more detail in the coming subsec- 
tions. Let us first assume that (j)(k;s,a) > for some 
positive k. For such values of k the normalized growth 
rate T reads 



r = \K] (A)|/A = |fc|V<Mfc;a,<7) 



(27) 



In the present paper we will consider situations where 
the MI of the nonlocal NLS equation is of finite band- 
width type and that each gain band has a unique maxi- 
mum growth rate, in the same as we have for the MI of 
the focusing NLS model and the weakly nonlocal limit, 
and it emerges as a generic property for a large class 
of response functions like the Gaussian, the exponen- 
tially decay function and the square pulse function. It is 
straightforward to derive the conditions which must be 
fulfilled: We find 



r' = \r 1/2 [20 + fc0'] 

and hence the maximum point k — fc max obeys 



(28) 



[20 + ^'] fc=fen 



= 



(29) 



The second derivative of T at k — fc max is given as 



T" = i<T 1/2 [30' + k$'\ 



k=k„ 



by N, is odd, i.e. N = 2m — 1 (to = 1,2, ..) and that 
there are to MI bands in this case. Notice that it al- 
ways exists a closed, bounded interval [0, ki] for which 
4>(ki) = and <j> (k) > for < k < k\. We refer to 
[0, fci] as the fundamental gain band. This means that 
we always have MI in the focusing case, due to the exis- 
tence of the fundamental band. In all the gain bands the 
growth rate is given by (p7j). 

Now, assume $ is a monotonically decreasing func- 
tion for all k > 0. In this case there is a unique 
fci > such that 0(fci) = 0, (j)(k) > (<j)(k) < 0) for 
< k < ki (fc > ki) i.e. there is only one band produc- 
ing MI. Now, since 0(0) = 1, <f>' (0) = and 0(fci) = 0, 
we find that T'given by |2|) satisfies T' (0) = 1 and 
lim fe _ >fe - T' (fc) = — oo. Hence, since </> and <fi' are con- 
tinuous, we can appeal to the intermediate value theo- 
rem and conclude that there is at least one k cr such that 

r' (k cr ) = o. 

It is also easy to study the variation of the growth rate 
curve with the width a. Simple computation yields 



d a T = hz 2 (t>- 1/2 & {ak) 



(33) 



and since $ is a monotonically decreasing function, it 
follows that the growth rate decreases with the width 
a for a fixed modulation wavenumber k. Moreover, the 
end point k\ of the MI band is a function of a, and the 
implicit function theorem yields the change of fci with 



(30) ^ 



dki 
da 



$' (ah) a -2k! 



(34) 



function <p must obey the inequality 



[30'+fc<n fc=fcn 



< 



(31) 



The condition ( pl| ) is referred to as the generic condition 
for the existence of a maximum growth rate. One can 
prove that if all the critical points of V obey (|3l|), these 
points coincide, and hence we have an MI pattern which 
resembles the MI of the local, focusing NLS. 



A. The focusing case (s=l) 

In this case the function defined by ( p2] ) appearing in 
(H) and (§?]) is given as 



4> = $ (ak) - k 2 



(32) 



From the list of general properties we find that <f> (0) = 1 
and we can prove that the number of zeros of <fi, denoted 



Since <&'(ak) < , we will have &(aki)a - 2fci < 
(which is the transversality condition </>' ^ evaluated 
at k = ki). Thus fci is a decreasing function of a, which 
means that the width of the MI band decreases in size. 
Let us use the Gaussian response function 



R(x) = U(-); $(C) = ^exp(~C 2 ) 
a \ a / y/ir 

as an example. The Fourier - transform is given by 

1 



R(k) = 3>(crfc) = exp 



- 2,2 

a k 

4 



(35) 



(36) 



which is a monotonically decreasing function. In Fig. 1 
it appears as a characteristic feature that the MI is be- 
ing suppressed as the characteristic widths a of these 
response functions are increased. Moreover, the MI gain 
band has a unique maximum growth rate. Notice that 
exactly the same characteristic features are found for the 
exponential respons<£3. 
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(5=1 



Modulational wavenumber k 
Fig. 1. The growth rate V as function of the wavenum- 
ber k for a — 1,2,5 and Gaussian response function in the 
focusing case. 

Now, let us consider situations where it is not assumed 
that the spectrum is strictly decreasing for all k. Then it 
is possible to have additional gain bands in the focusing 
case. We assume that 



<j>{k 2 j) = 4>{k 2o +i) = 0; 



(f> (k) > for k e (k 2 j, k 2 j+\) 



(37) 



(38) 



for j — 1, 2, ..,771— 1 where m is the number of MI 
bands. Since T = kcj) 1 ^ 2 is continuous and realval- 
ued on [k 2 j, fcaj+i] an d differentiable on (k 2 j,k 2 j+i) , it 
follows from Rolles theorem that there is at least one 



fccr £ (k 2 j, k 2j+ i) such that T'(kci') = 0. If, in addition 



(Jh 



<f>'{k 2j ) > 0, tffaj+i) < 0, 
we find the following limits 



(39) 



lim r (k) = +oo, 

k— >k 2 j 



lim 



r'(fc) 



Now, if the condition ( |3l| ) is satisfied, then k. 



Si) 



is a unique. Hence we have a unique maximum growth 
rate for this wavenumber band as well. 

As an example we consider the square pulse function 



i?(x) = !$(-); 0>(C) 



h Kl<i 

0; |CI>1 



whose Fourier - transform is given by 
sin(fccr) 



iZ(ft) = $(crfc) 



ka 



(40) 



(41) 



For small and moderate values of the width a we only 
have one fundamental MI gain band. The bifurcation 
equation reads 

(42) 



tan(z) = z/3 



in this case, and by means of fl25p with $(z) > 
sin(z) > 0), we find the critical values a c . More- 
over, one finds that <\>" — — z 2 — 6, from which it follows 
that new MI bands are excited at the bifurcation points. 
The table below summarizes the findings for the lowest 
order excitations. 



a c ,i = 21.03 


The second MI band is excited 


er C) 2 = 52.555 


The third MI band is excited 



(43) 



A simple computation based on ( |3i"| ) now reveals that 
each MI band has a maximum growth rate. Finally, 
since $ of the square response function is a monotoni- 
cally decreasing function on an interval containing 0, hi] 
for all values of a, we conclude by appealing to (|33) and 
( |34] ) that an increase in the width a decreases both the 
fundamental band and the maximum growth rate. The 
expressio n (p4| ) (with k\ replaced with ki, (i = 2j, 2j+l)) 
and (|37|)-(|39|) show that the higher order MI gain bands 
move towards lower wavenumbers as the width a in- 
creases. Fig. 2 summarizes these features graphically, 
and the results are consistent with the findings ii£3. 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

Modulational wavenumber k 
Fig. 2. The growth rate T as function of the wavenumber 
k for a — 20, 40, 90 and the square response function in the 
focusing case. 



B. The defocusing case (s=-l) 

The function (f> given by 
4> = -%(ak) - k 2 ; 



(44) 
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in the defocusing case. Here the number of zeros of <f>, 
denoted by N, is even, i.e. TV = 2to (m = 0, 1, 2, ..) and 
there are to MI bands. We observe from the list of prop- 
erties of the spectrum $ (ak) that MI can only exist in 
the high wavenumber regime, as opposed to the focus- 
ing case. Notice that the case to = corresponds to 
the situation where <f) (k) < for all k > which means 
modulational stability. 

A typical situation revealing stability occurs when 
$(crfc) > for all k > 0, irrespective of the value of a. 
The spectrum of the Gaussian ([36|) satisfies this prop- 
erty. The same holds true for the exponential response 
function 



R(x) = U(-) ] 3(0 = iexp(-|C|) 
a \<jJ 2 

for which the Fourier - transform is given by 

1 



R(k) = <f>(ak) = 



1 + <7 2 fc 2 ' 



(45) 



(46) 



When $(fc) changes sign, the picture is somewhat 
more complicated. First of all, we will still have stability 
if $(erfc) < for all k provided the control parameter a is 
below certain threshold a c . The threshold value a c and 
the corresponding value of k can be found by requiring 
the bifurcation condition (^3|) to be fulfilled. Secondly, if 
<j belongs to the complementary regime, i.e., a > a c we 
have <j>(crfc) > for certain wavenumber intervals, and 
hence we may have one or more MI gain bands, whose 
growth rate T is given by (^) and ([44]). Even in this 
case one can prove the existence of a unique maximum 
growth rate within each unstable band in the same way 
as in the focusing case. 

Again we use the square pulse function ( p7j| ) - ([ll]) as 
an example. The bifurcation problem is solved by means 
of ([l2]) - (^B|) with s = — 1 and hence the extra constraint 
that sin (z) < 0. Here we also find that <j>" = — z 2 — 6 < 0, 
which means that new MI bands are excited through the 
bifurcation process. The table below displays the bifur- 
cation values of a for the lowest order excitations. 



cr c ,i = 9.146 


The first MI band is excited 


cr c ,2 = 35.78 


The second MI band is excited 


<Jc,3 = 71.31 


The third MI band isexcited 



(47) 



The number of MI bands increases when the width in- 
creases. By ([H]) each MI band has a maximum growth 
rate. Finally, by using the expression for the velocity 
of the zeros fej-i and k^j of <f> (j = 1,2,.., to) in the 
defocusing case 



dki &(aki)ki 
da (aki) a - 2ki 



2j-l,2j 



(48) 



and (|37|) - (£59|) with <j> given as (J4J) and k 2 j and 
^2j+i replaced with A^-iand ^ , respectively, it can be 
shown in a the same way as in the focusing case that 
^r 1 < and hence the MI bands move to the low 

da 

wavenumber regimes as the width increases. The latter 
analytical result is also,in accordance with the numeri- 
cal results obtained irE3. In Fig. 3. the MI result for the 
square pulse (M) in the defocusing case is displayed. 



s 

o 

(5 04 




a = 20 

~ - o = 40 
o = 90 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

Modulational wavenumber k 
Fig. 3. The growth rate F as function of the wavenum- 
ber k for a — 20, 40, 90 for the square response function in 
the defocusing case. 



3. Generalized response functions 

Now, let us consider a localized response function R(x) 
which can be written as a convex combination of N re- 
sponse functions R n (x) whose properties are known i.e. 
the weighted mean 



A I 



(49) 



71 = 1 



where 



M 



/oo 
R n (x)dx = l; 

(50) 

and R n {x) is parametrized by a typical lengthscale a n , 
i.e. R n (x) = a^^n with as the scaling func- 

tions. From now on we refer to the functions R n (x) as 
" building blocks" . 

The corresponding spectrum R (k) appears as a con- 
vex combination of R n (k) 
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M 



(51) 



where $„ is the spectrum of the scaling function $ r , 
Moreover, we find in accordance with (H) that 



x 2 R n (x) dx ~ rj^; 



(52) 



A typical lengthscale <r of the composite response func- 
tion R (x) can be defined 



iV 



^ 2 = ^ n<7 ™ ; 



(53) 



n=l 



It is straightforward to show that the dispersion relation 
on normalized form can be expressed as 



where 



A I 



s 5Z ^ n$n ( ank ) _ fc2 ' 



(54) 



(55) 



which is a generalization of the dispersion relation (JlS)) . 
Here Ao and fe are given as ( p0|) and (21), respectively, 
while cr„ is the redefined width of the response function: 



(56) 



P0&n — * 0"n 



Now, if the transversality condition 4>' conv ^ is satisfied 
on points where 4> C onv = 0, there are a finite union of 
disjoint, bounded and closed subintervals of the positive 
k axis, for which <\> conv > 0. This means that the MI is of 
finite bandwidth type in this case as well. In each gain 
band the normalized growth rate T is given as 



|ftl(A)|/Ao = |fc|^ 



(57) 



In (54) - ( p7|) there are 2M — 1 control parameters, 
namely the M width parameters g\,<j%, ...,<tm and the 
M — 1 weight parameters fi , fa , . . . , Jm- i ■ To summarize, 
the stability properties depend on the weight /„ of each 
Fourier - transform <!>„ and length scale o~ n . 

Let us consider some special cases. Assume that the 
Fourier transforms R n (k) of all the building block func- 
tions R n {x) are monotonically decreasing positive func- 
tions, such as the exponential and Gaussian response 
functions. The sum J2 n =i fn&n(°~nkt) will now be a de- 
creasing and positive function of k. In the defocusing 
case (s = — 1) we will have modulational stability in this 
case, while we get MI with one gain band about k = in 
the focusing case (s=l), which has a unique maximum 



growth rate. In Fig. 4. we display this phenomenon for 
the case when the total response function is a sum of a 
Gaussian and an exponential response function. 




Modulational wavenumber k 
Fig. 4. Growth rate F versus wavenumber k for a con- 
vex combination of a Gaussian response function with weight 
factor f g = / and width u g — l and an exponential response 
function with weight factor f e — 1 — f and width a e =2 for 
different values of /. 



In Fig. 5 we show what happens for a focusing non- 
linearity when the response function is a 3-component 
combination of a Gaussian, an exponential, and a square 
response function. In Fig. 5 two MI gain bands are re- 
vealed. A typical feature of the fundamental band is that 
the growth rate exhibits a more complex behavior with 
several local extrema, as opposed to the single-maximum 
generic situation defined by (|31 




0.15 0.2 0.25 

Modulational wavenumber k 



7 



Fig. 5. Growth rate P versus wavenumber k for a con- 
vex combination of a Gaussian (width a g =10), an exponen- 
tial (width a e —20), and a square response function (width 
a sq =90) in the focusing case. The corresponding weight fac- 
tors are f g —0.1, / e =0.3, and / S9 =0.6. 

4. Conclusion 

The linear stage of the MI for the nonlocal NLS equation 
is conveniently studied in terms of the spectrum of the 
response function. From the dispersion relation ([l7]) it 
appears that the crucial point in this discussion is the 
location of the spectrum of the response function rela- 
tive to the parabola in fc-space. The following features 
complement and extend the results obtained inB3: 

• The MI is of the finite bandwidth type. It consists 
of a finite number of well-separated gain bands. 
Moreover, it is possible to predict the occurrence 
of excitation, vanishing, coalescence and separa- 
tion of MI bands. 

• For a large class of response functions, each MI 
band has a unique maximum growth rate. This 
property holds true for the Gaussian, the exponen- 
tial and the square response functions. It resembles 
the structure of the MI bands in the focusing local 
NLS equation. 

• In the focusing case we always find at least one 
MI gain band centered about k = 0. It is ver- 
ified analytically that the width of this MI band, 
as well as the corresponding growth rate, decreases 
when increasing the width of the response function, 
provided the spectrum of the response function is 
decreasing in this MI band. Furthermore, addi- 
tional MI bands are excited at higher wavenum- 
bers when the width parameter exceeds a certain 
threshold, i.e. when the nonlinearity becomes suf- 
ficiently nonlocal. The latter phenomenon is a 
unique feature of the nonlocal nonlinearity and has 
no equivalent in the local case and the weakly non- 
local limit. 

• In the defocusing case we can either have stabil- 
ity or MI of the finite bandwidth type. The latter 
situation can only occur in the high wavenumber 
regime, and only if the width of the response func- 
tion exceeds a certain threshold, i.e. when the non- 
linearity becomes sufficiently nonlocal. 

• In both the focusing and defocusing case the higher 
order MI bands move towards lower wavenumbers 
as the width of the response function increases. In 



the limit of strong nonlocality the MI bands van- 
ish completely. This result agrees with the fact the 
strongly nonlocal limit of the NLS model ([!]) is a 
linear model. 

• It is possible to study MI for more complicated 
multiscale scenarios where the response functions 
decompose into a weighted mean of "building 
blocks" , where each " building block" has a char- 
acterized width. In this case we get a richer reper- 
tory of response functions to deal with. This aims 
at showing that the results we have obtained rep- 
resents generic features of the MI of the nonlo- 
cal NLS with general symmetric, positive response 
functions. 
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